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ABSTRACT 


Tne  advent  of  jet  fighter/attack  aircraft  has 
resulted  in  the  need  for  additional  devices  to  aid  the 
pilot  in  the  accurate  delivery  of  unguided  weapons. 
The  current  generation  of  attack  aircraft  has  met  this 
need  with  the  installation  of  minicomputers  and 
additional  sensor  devices  such  as  doppler  radar  and 
inertial  navigation  systems. 

This  paper  demonstrates  a  low  cost  method  for 
improving  the  delivery  accuracy  of  those  aircraft 
which  do  not  have  such  devices  installed.  A 
comparison  is  made  between  delivery  accuracy  with  and 
without  the  proposed  modification. 
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I.   INTRODUCTION 


The  purpose  of  this  paper  is  to  examine  the  feasibility 
of  using  a  micro  computer  in  Naval  aircraft  to  aid  the  pilot 
in  the  solution  of  the  ballistics  problem  for  unguided 
weapons.  Specifically,  this  problem  is  to  release  a  weapon 
at  some  point  in  space  with  the  proper  velocity  vector  so 
that  the  weapon  will  eventually  strike  a  selected  target. 
Since  several  models  of  attack  aircraft  (A-6E,  A-7E)  are 
equipped  with  computers  which  perform  this  function,  special 
attention  will  be  focused  oa  those  aircraft  in  which  the 
pilot  manually  provides  the  solution.  Such  aircraft  are 
referred  to  as  fixed  sight  aircraft  and  the  F-4  is 
representative  of  this  class. 

The  micro  computer,  because  of  its  low  cost  and  weight, 
is  a  device  with  the  potential  to  provide  the  pilot  of  fixed 
sight  .aircraft  with  useful  and  timely  information.  To 
determine  whether  such  a  potential  could  be  exploited,  the 
INT3L  8080  is  examined  as  being  representative  of  modern 
micro  computers  and  the  F-^4  is  considered  as  a  typical  fixed 
sight  aircraft.  The  concepts,  however,  are  intended  to  be 
more  general  in  nature  and  no  discussion  ,  data,  or 
conclusion  contained  herein  is  intended  as  an  endorsement  or 
criticism  of  the  manufacturer  or  weapon  systen  involved. 

The  body  of  the  paper  consists  of  four  major  sections. 
Tiie  sections  are  as  follows:  a  description  of  the  general 
nature  of  the  unguided  weapons  delivery  problem  and  how  the 
problem  is  now  solved  by  pilots  of  fixed  sight  aircraft;  the 
proposed   use  of   a   micro   computer   to   aid  the  pilot;  an 


analysis  and  test  of  the  proposed  solution;  and  conclusions 


II. 


ONGUIDED     WEAPONS    DELIVERY    USING    FIXED    SIGHT    AIRCRAFT 


A. 


DEFINITIONS 


Fixed  sight  aircraft  is  a  term  applied  tc  these  aircraft 
equipped  with  neither  a  lead  computing  gunsight  nor  a 
computer  operated  bomb  release  device.  In  such  aircraft  a 
sighting  device  or  gunsight  is  located  in  the  pilot's 
forward  field  of  view.  Its  purpose  is  to  represent  visually 
an  angle  in  milliradians  (mils)  measured  in  the  positive 
direction  'below  a  longitudinal  reference  line  called  the 
FUSELAGE  REFERENCE  LINE  (FRL)  .  Figure  1  depicts  such  a 
device.  Typically  these  devices  are  adjustable  through  a 
range  of  zero  to  approximately  200  mils,  but  once  adjusted 
they  remain  set  until  manually  changed  by  the  pilot. 


PARALLEL  BORESIGHTED  I  un FUSELage refersSI^ 

L  FUSELAGE  ANGLE  OF  ATTACK 


Fig  1   SIDE  VIEW  AND  PILOT  VIEW  OF  A  GUNSIGHT 
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B»    P30BLEH  DEFINITION 


For  a  particle  of  constant  mass,  m,  moving  under  the 
action  of  an  applied  force,  F,  Newton's  second  law  of  motion 

states  that 

F=m*a.        (1) 

If  the  weapon  is  considered  to  be  such  a  particle,  and 
if  the  position  of  the  weapon  is  measured  in  a  rectangular 
coordinate  system,  (x,y,z),  then  equation  (1)  may  be  written 
in  vector  form  as 


(a  ,a  ,  a  )=   (F  ,F  ,F  )/m.   (2) 
x   y    z      x   y   z 


For   air   to   ground  weapons   delivery,   the  coordinate 

system  is  chosen  so  that  y  represents  altitude,  x  represents 

the  down  range  travel  of  the  weapon,  and  the  z  axis  measures 

displacement  in  azimuth.  Ihe  origin  is  chosen  so   that,   at 

the  moment  of  release,  x  equals  zero. 

If  it  is  assumed  that  the  air   mass   through   which   the 

weapon   is   moving   is  stationary  with  respect  to  the  earth, 

that  gravitational  forces  are  constant  with   respect   to   y, 

and   that   forces   induced   by  the  rotation  of  the  earth  are 

negligible,  then  F  =0.   The  remaining  forces  acting   on   the 

z 

weapon   are   drag,  thrust,  and  gravitational  force.   It  will 

be  assumed  that  thrust  is  zero,  for  all  weapons  except 
rockets.  Rockets  are  assumed  to  have  a  short  initial  period 
of  constant  thrust  after  which  they  fall  with  the  same 
discipline   as  all  other  unguided  weapons.   Jrag  acts  in  the 
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direction  opposite  to  the  direction  of  the   velocity   vector 

2 
and   is  directly  proportional  to  ||V||   and  air  density.   In 


a  standard  air  mass,  density  is   a   decreasing   function   of 

altitude.     Therefore    equation   (2)   may   be   written   in 
component  form  as  follows 


a  =-drag  *  1/m     (3) 
x      x 


a  =   -(drag  +g*m)/m   (4) 
y        Y 


2 

where  g  is  a  constant  32.17  ft/sec  . 


Letting  THETA=ARCTAN (v  /v  ),  (3)  and  (4)  become 

y  x 


and 


a  =-drag*cos (THETA) /m     (5) 
x 


a  =-drag*sin  (THSTA) /ra   -g.    (6) 

y 


Drag  is  represented  as 


2 

drag=rho*S*V  *C 

d 

where  'rho'  is  the  air  density,  'S'  is  the   cross   sectional 

2    2  1/2 
area,   ' V   is   the  total  velocity  (  (V   +V   )     ) ,  and  C  , 

x    y  d 

the  drag  coefficient,  is  a  dimensionless  factor  whose   value 

depends  on  the  shape  of  the  weapon  and  the  ratio  of  *V'  to 
the  speed  of  sound  at  altitude  y.  Equations  (5)  and  (6)  may 
be  considered  as  differential  equations 
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a  =f(v  ,v  ,x,y,t)     (7)  and 
x     x   y 


a  =h(v  ,v  ,x,y,t)  .    (8) 
Y     x   y 


In  air-to-ground  weapons  delivery,  equations  (7)  and  (3) 

may   be  solved  numerically  for  any  set  of  initial  conditions 

0   0   0   0 
(v  ,v  ,x  ,y  )  to  determine  a  time  t1   where   y   =the   target 
x   y  t « 

elevation   and   x     is  the  down  range  travel  of  the  weapon. 

t« 

Therefore,  the  problem   for   air-to-ground   delivery   is   to 

position   the   aircraft   so   that  the  initial  conditions  are 

matched  when  the  selected  aim  point  is  a  distance   x     down 

t« 

range  from  the  aircraft. 


For  unguided  air-to-air  weapons,   equations   (7)   and (8) 
may   be   solved   together   with   initial  value  conditions  to 
describe  the  position  (x  (t)  , y  (t)  , z  (t) )  of  the  weapon  at   any 
time   t  after  the  weapon  has  been  fired.   Since  the  position 
of    the    target    is    also    some    function    of    time, 

(x       (t),y       (t),z       (t)  )  ,  the  problem  is  to  find  a 
target     target     target 

set  of  initial  value  conditions  such   tnat   for   some   time, 

t'«, 

U(t")  ,y  (t")  ,z  (t"))  =  (x       (f),y       (t««),z       ( 

target       target       target 

t"))  . 


In  either  case,  the   aircraft   instruments   provide   the 


13 


indication  that  the  initial  conditions  are  satisfied.  The 
gunsight  is  used  to  establish  a  spatial  relationship  to  the 
target.  The  desired  terminal  condition  occurs  at  so;ae 
future  time  when  the  weapon  hits  the  target.  This  mixture 
of  initial  condition  and  terminal  condition  constraints 
makes  the  problem  a  combination  initial  value  problem  and 
boundary  value  problem. 


C.    PILOT  SOLUTION  FOR    AIR  TO  GROUND  DELIVERY 


Although  high  speed  computers  which  numerically  solve 
(just  prior  to  weapon  release)  equations  (7)  and  (8)  are 
installed  on  some  attack  aircraft  (A-6E,  A-7E) ,  no  such 
devices  are  installed  on  fixed  sight  aircraft.  For  the 
pilot  of  the  fixed  sight  aircraft  a  numerical  solution  of 
the  problem  is  not  available  in  flight.  Therefore  sets  of 
solutions  have  been  calculated  for  various  initial 
conditions.  The  pilot  decides,  prior  to  reaching  the 
target,  which  of  the  precomputed  sets  of  initial  conditions 
he  will  match.  Then  he  consults  a  table  which  contains  a 
standard  sight  setting.  After  adjusting  this  setting  to 
correct  for  variations  in  aircraft  gross  weight  and  release 
altitude,  he  adjusts  his  sight  to  this  corrected  value. 
Finally,  he  flies  the  aircraft  so  as  to  match  the  delivery 
parameters,  places  the  sight  image  over  the  target,  and 
releases  the  weapon.  If  all  parameters  have  been  matched  , 
the  weapon  will  have  the  greatest  probability  of  striking 
the  target. 

In  actual  deliveries,  however,  the  parameters  are 
frequently  not  matched  exactly.  This  is  understandable 
considering  that  the  controls  available  to  the  pilot 
(throttle  and  elevator  control)  do  not  directly, 
instantaneously  or  independently   control   the   position   or 
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initial  velocity  vector.  Furthermore,  dive  angle 
information  is  not  directly  available  from  any  aircraft 
instrument  and  the  information  which  is  directly  available 
(true  airspeed  and  altitude)  must  be  mentally  processed  by 
the  pilot  before  he  can  physically  react  ,  thereby  adding 
further  delay.  Finally,  the  aircraft  itself  cannot  react 
instantaneously  to  control  inputs.  Gauthier^/"-!^  has  shown 
that  sudden  perturbations  of  the  angle  of  attack  caused  by 
pilot  induced  changes  in  pitch  attitude  (  elevator  control) 
will  cause  the  information  provided  by  a  fixed  gunsight  to 
be  unreliable  for  several  seconds.  In  addition  ,  as  shown 
in  Fig  2,  the  gunsight  image  is  constantly  moving  across 
the  earth  so  that  the  proper  sight  'picture1  occurs 
precisely  once  during  the  dive. 
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8000' 


14000' 


NOTE:     Although  Depression  from  flight  path    (106  mils )   extracted 
from  the  ballistic   tables  is  corrected  for  parallax,    the 
remaining  calculations   ignore   these  factors.^     There jore 
-it  is  apparent   that  there  are  small  geometric  anomalies 
in  the  above  diagram. . 


Fig    2       SIGHT    LINE    MOTION    IN    A    DIVE 
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In  spite  of  these  difficulties,  there  are  strategies 
which  help  the  pilot  to  learn  a  solution  to  the  problem. 
First,  the  deliveries  have  been  standardized  so  that  the 
pilot  may  concentrate  his  efforts  at  mastering  several  basic 
maneuvers.  This  practice  tends  to  reduce  the  magnitude  of 
delivery  errors.  The  second  aspect  which  aids  the  pilot  is 
known  as  error  analysis.  Since  actual  conditions  rarely 
match  those  cf  the  practice  patterns,  the  pilot  is  taught  to 
recognize  those  errors  which  are  present  and  to  apply 
corrective  measures.  Fig  3  graphically  depicts  the  effect 
of  various  errors  in  airspeed,  dive  angle,  release  altitude, 
and  aim  point.  For  each  of  the  standard  deliveries,  there 
is  a  table  in  reference  / 2 _f  which  lists,  for  each  weapon, 
the  effect  upon  weapon  impact  of  errors  of  10  knots  in  true 
airspeed  (  fast  arid  slow  )  ,  one  degree  in  dive  angle  (  steep 
and  shallow),  100  feet  in  release  altitude  (  high  and  low), 
and  one  raiilirad  ian  in  sight  setting  (  plus  and  minus)  . 

To  use  the  tables,  the  pilot  memorizes  key  values  for 
the  weapons  he  is  going  to  release.  During  the  tracking 
portion  of  his  dive,  he  first  attempts  to  minimize  the 
errors  in  each  of  the  delivery  parameters.  If  he  decides 
that  the  value  of  some  or  all  of  the  delivery  parameters 
will  not  egual  the  preselected  Values,  he  estimates  the 
cumulative  effect  of  all  such  errors  on  weapon  impact. 
Having  estimated  this  miss  distance,  h<  then  deliberately 
causes  another  error  either  in  release  i ltitude  or  sight 
angle  to  counteract  the  effect  of  the  previous  errors. 
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SLOW 


LONG 


SHORT 


RELEASE  HEIGHT  ERROR 


LEVEL    DELIVERY 

HIGH  RELEASE   HEIGHT 


PLANNED  RELEASE   HEIGHT 


*.  LOW  RELEASE   HEIGHT 


SHORT  LONG 


DIVE  ANGLE  ERROR 

CONSTANT   RELEASE    HEIGHT  AND  AIRSPEED 


SHORT 

SIGHT  ANGLE  ERROR 


LONG 


LONG 


SHORT 


Fig  3   EFFECTS  OF  AIRSPEED,  HEIGHT,  DIVE,  AND  SIGHT  ERRORS 
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As  an  illustration  of  this  method,  suppose  the  pilot 
estimates  that  he  will  be  10  knots  fast  and  one  degree  steep 
at  release  altitude  and  that  the  combined  effect  of  the 
errors  will  place  the  weapon  100  feet  past  the  target.  If, 
in  the  error  sensitivities  chart  for  the  particular  weapon, 
the  entries  for  +100  ft.  and  -1  mil  are  -25  feet  and  -10 
feet  respectively,  then  the  pilot  may  choose  to  release  the 
weapon  400  feet  high  in  relation  to  the  preselected  release 
altitude  or  he  may  use  an  imaginary  sight  line  which 
represents  an  angle  of  10  mils  less  than  the  preselected 
angle. 

In  view  of  the  difficulties  involved  in  fixed  sight  dive 
bombing,  it  is  not  surprising  that  considerable  effort  has 
been  devoted  to  the  development  and  implementation  of 
computer  systems  to  improve  the  delivery  accuracy  of  all 
current  generation  N"aval  Attack  aircraft. 


D.    AIR  TO  AIR  GUNNERY 


In  air-to-air  gunnery,  the  pilot  of  a  fixed  sight 
aircraft  essentially  solves  his  aiming  problem  in  the  same 
manner  as  for  air  to  ground  weapons  delivery.  He  selects 
several  representative  situations  and  consults  a  table  to 
determine  the  correct  sight  setting  for  each  particular 
situation.  In  air  to  air  gunnery,  however,  the  pilot  would 
usually  like  to  be  able  to  fire  his  aircraft's  guns  at  any 
time  he  has  a  good  chance  of  hitting  the  target.  Therefore, 
rather  than  try  to  match  any  preselected  set  of  firing 
conditions,  the  pilot  takes  the  best  guess  and  fires.  The 
guns  are  loaded  with  a  mixture  of  tracer  rounds  so  that  the 
pilot  visually  adjusts  his  aim  depending  on  his  evaluation 
of  the  effect  of  previous  aim  angles. 
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III. 


PROPOSED  SOLUTION 


In  order  to  achieve  the  full  benefit  of  the  low  cost  of 
the  micro  computer,  it  was  assumed  that  no  additional  sensor 
equipment  such  as  doppler,  air-to-ground  ranging  devices,  or 
inertial  navigation  devices  would  be  included  in  the 
modified  system.  Instead,  the  system  would  be  designed 
around  existing  components  with  ultimate  responsibility  for 
problem  solution  still  resting  with  the  pilot.  The  micro 
computer  is  used  chiefly  to  perform  those  calculations  which 
the  pilot  of  fixed  sight  aircraft  now  estimates. 


A. 


COMPONENTS 


The  components  of  the  modified  system  are  depicted  in 
Fig  4.  Their  relationship  to  the  overall  system  is  given 
below. 
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1 .    Control  S tick  Sense  Switch 

An  on/off  switch  which  is  activatad  'on'  when  the 
pilot  begins  tracking  the  target. 

2 •   £o£j;isble  Micro  Corapu ter  With  A  uxil  iary_  Mem  o  ry_ 

A  device  whose  auxiliary  storage  contains  a  complete 
library  of  standard  release  parameters  and  error  sensitivity 
factors  for  each  weapon  authorized  for  carriage  and  release 
from  the  particular  aircraft.  Each  weapon  would  require 
approximately  128  bytes  of  storage.  It  is  assumed  that  no 
more  than  50  such  weapons  would  apply  to  any  particular 
aircraft.  The  ground  based  micro  computer  would  be  capable 
of  retrieving  blocks  of  such  data  and  passing  the  data  to 
the  micro  computer  aboard  the  aircraft.  It  is  also  assumed 
that  no  more  than  three  different  weapons  would  be  carried 
by  a  particular  aircraft  on  any  flight. 

3«    Air  Ba.ta  Computer 

A  device  currently  installed  on  Naval  aircraft  which 
provides  analog  signals  corresponding  to  true  airspeed, 
calibrated  airspeed,  and  altitude.  These  signals  would  be 
converted  to  digital  format  to  be  read  by  the  micro 
computer. 

**  •   Loft  Bombi  ng  Computer 

An  attitude  reference  gyro  of  the  type  currently 
installed  in  the  F-4  whose  analog  signals   corresponding   to 
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pitch  and  roll  angle  as  well  as  rate  of  pitch  and  heading 
change  would  be  converted  to  digital  format  to  be  read  by 
the  micro  computer. 

5.  Cockpit  Keyboard 

Any  standard  numerical  keyboard,  of  the  type  used  on 
pocket  calculators,  which  has  at  least  five  function 
switches  and  light  emitting  diodes  to  confirm  register 
contents.  This  device  would  be  used  by  the  pilot,  or 
preferrably  by  a  crewman,  to  input  the  mode  (  air-to-air  or 
air-to-ground) ,  and  if  in  the  air-to-ground  mode  to  input 
the  desired  delivery,  weapon  to  be  released,  target 
elevation,  and  aircraft  gross  weight. 

6.  Micro  Computer 

A  micro  computer  with  a  word  size  of  at  least  eight 
bits  and  speed  at  least  that  of  the  INTEL  8080  (2.5  micro 
second  instruction  cycle) . 

7.  Gun sight 

Either  a  standard  fixed  gunsight  modified  by 
replacing  the  fixed  reticle  light  source  with  an  array  of 
addressable  light  emitting  diodes,  or  a  standard  'heads  up' 
gunsight. 

8«   Airborne  Missile  Control  System 

One  of  the  several  models  of  air-to-air  radars 
currently  installed  on  the  F-4  which  provides  analog  signals 
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corresponding  to  the  present  range  to  the  target  and  the 
rate  of  closure.  These  analog  signals  would  be  converted  to 
digital  format  to  be  read  by  the  micro  computer. 

9.   AH<ll~  21.  kll&ck    Sensor 

A   device   which   is  currently  installed  on  Navy  jet 

aircraft  which  provides  an  analog   signal   corresponding  to 

the   angle   between   the  flight  path  of  the  aircraft  and  the 

FUSELAGE  REFERENCE  LINE.  This  signal  would  be  converted  to 
digital  format  for  use  in  air-to-air  calculations. 


B. 


DESCRIPTION  0?  USE 


Fig  5   shows   the   overall   timing   sequence   for    the 
functions  of  the  system. 
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1  •   Air  to  Ground  ^ode 

Prior  to  flight,  the  aircraft  is  loaded  with  some 
authorized  ordnance  configuration.  After  the  aircraft  is 
operating  on  internal  electrical  power,  a  portable  micro 
computer  loads  the  appropriate  tables  into  the  memory  of  the 
aircraft's  computer.  As  the  aircraft  approaches  the  target, 
the  pilot/crewman  keys  a  standard  delivery  and  weapon  code 
into  the  micro  computer  along  with  the  target  elevation  and 
aircraft  gross  weight.  As  the  pilot  commences  his  roll-in, 
the  micro  computer  samples  calibrated  airspeed  and  pitch 
angle  and  calculates  the  angle  of  attack  of  the  fuselage 
reference  line  for  a  steady  state  dive  at  that  pitch  angle 
and  calibrated  airspeed  and  displays  the  angle  of  attack  in 
mils  on  the  gunsight.  The  pilot  uses  this  information  to 
establish  the  aircraft  in  a  dive  patn  which  intersects  the 
earth  at  some  predetermined  distance  beyond  the  target  (  aim 
off  point) .  When  the  aircraft  is  established  in  the  proper 
dive,  the  pilot  actuates  the  control  stick  sense  switch  and 
the  coaputer  begins  to  sample  true  airspeedr  calibrated 
airspeed,  altitude,  and  pitch  angle.  The  pilot  makes  no 
farther  throttle  adjustments.  The  computer  calculates  dive 
angle  and  estimates  what  the  release  parameters  {  airspeed 
and  dive  angle)  will  be  when  the  aircraft  reaches  the 
prescribed  release  height  above  the  earth.  This  sampling 
and  calculation  takes  about  two  seconds.  The  computer  then 
calculates  the  effect  of  any  errors  in  the  estimated 
airspeed  and  dive  angle  and  corrects  for  the  impact  error  by 
addressing  a  new  diode  in  the  gunsight  or  sending  a 
correction  signal  to  the  heads  up  gunsight.  The  new  sight 
setting  would  be  calculated  using  the  error  sensitivity 
factors  which  had  been  read  prior  to  flight.  The  pilot  uses 
the  corrected  sight  setting  as  the  new  aiming  reference  and 
places  the  revised  sight  pipper  over  the  target   at   release 
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altitude  arid  manually  releases  the  weapon.  After  dive 
recovery,  the  computer  is  reset  and  the  process  is  repeated 
as  desired. 

2 •   Air  12  Air  Mode 

For  air  to  air  use,  there  would  be  no  need  to  load 
tables  into  the  aircraft  computer  prior  to  flight.  Instead, 
the  program  and  appropriate  tables  would  remain  in  the 
computer  memory. 

In  flight  the  pilot/crewman  would  select  the 
air-to-air  mode.  The  computer  would  then  begin  to  sample 
true  airspeed,  altitude,  present  range  to  the  target, 
closure  rate,  roll  angle  and  pitch  and  azimuth  rates  of 
change. 

At  this  time  the  35  mil  setting  would  be  illuminated 
on  the  gunsight.  w'hen  the  pilot  had  established  a  suitable 
tracking  position  with  the  sight  held  steady  on  the  target, 
he  would  activate  the  control  stick  sense  switch. 

The  computer  would  then  calculate  the  proper  lead 
angle  for  the  most  recently  sampled  conditions  and  display 
this  lead  angle  on  the  gunsight.  The  pilot  would  adjust  to 
the  new  aim  angle  as  smoothly  and  quickly  as  possible  and 
fire  the  guns.  Tracers  would  still  be  mixed  with  other 
ammunition  so  that  the  pilot  could  adjust  his  fire  as 
necessary.  The  time  required  to  perform  the  calculations 
would  be  about  .2  seconds. 
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IV.    ANALYSIS  AND  TEST 


A.   AIR  TC  GROUND 


1 .   General  Consi deration s 

If  it  is  assumed  that  density,  temperature,  and  the 
speed  of  sound  are  continuous,  diff erentiable  functions  of 
altitude,  then  there  exist  unique  solutions  of  equations  (7) 
and  (8)  for  each  set  of  initial  conditions.  Moreover,  it  is 
intuitively  clear  that  if  the  solutions  were  written  in  the 
component  form 


0   0   0   0 
x(t)=F(v  ,v  ,x  ,y  ,t)   (9) 

x   y 

and 


0   0   0   0 
y(t)=H(v  ,v  ,x  ,y  ,t)  ,  (10) 
x   y 


then   the   position   of   the   weapon   at   time   t   could   be 

considered   to   be   a   continuous   function   of   the  initial 

conditions.   That  is,  for  any  neighborhood,  N,   of   a   point 

(x  (t) ,Y  (t) )   which   represents  the  position  of  the  weapon  at 

time,    t,    after    release    (with    initial    conditions 
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0   0   0   0 
(v  ,v  ,x  ,  v  )),   there   exists   a   neighborhood,   J,   of  the 
x   y 

initial  conditions  such  that  any  identical   weapon  released 

with   initial   conditions  which  were  in  j"  would,  at  tire,  z, 

be  in  H.   Furthermore,  it  aay  be  assumed  that  first  partial 

derivatives  of  F  and  H  exist  within  D  so  that 


0  0  0  0 

dX=F    *   dV    +   F   *  dV    +   F   *  dX   +   F  *  dT   (11) 

v       x      v      y      x  y 

x  y 


and 


0  0  0  0 

dY=H    *   dV  +   H    *  dV    *   H   *  dX   ♦   H   *  dY.  (12) 

v        x       v       v       x  y 

x  y 

Where  the  variables  and  partials  en  the  right  represent 
changes  in  the  initial  conditions  and  those  on  the  left 
represent  position  at  some  time,  t,  which  is  held  constant  . 

Since  the  proposed  system  incorporates  no  additional 
sensor  devices,  it  Gust  perform  essentially  the  same  type  of 
calculations  which  the  pilot  of  a  fixed  sight   aircraft   now 
performs.   It  samples  airspeed,  altitude,  and  pitch  attitude 
and  predicts  what  the  value  of  true  airspeed  and  dive   a:.  - 
will  be  at  release  altitude.   These  values  are  then  ce 
with,  the  standard   release   parameters   for   the   particular 
delivery   being   used.   Any  difference  between  the  predicted 
value  and  the  desired   value   is   converted   to   a   gunsight 
correction   which  must  be  applied  to  the  sight  to  counteract 
the  effect  of  the  release  airspeed  or  dive  angle  error, 
sum   of   the   correction  for  airspeed  and  the  correction  for 
dive  angle  is  then  applied  to  the   gunsight   and   the   pilot 
flies   the   aircraft   to   place  the  new  sight  image  over  the 
target  at  release  altitude. 

Since   the   proposed   solution   relies  on    the  error 
sensitivity  tables,  ana  since  these  tables  are  based  on   the 


partial   derivative   of   F   with   respect   to  .several  of  the 

components  of  the  initial  condition  {  for  example,  an   error 

in   true   airspeed   changes   v    ,v  ,and   x  ),  and  since  the 

x     y 

solution  assumes  that  effects  of  changes  in   airspeed,   dive 

angle,  and  aim  angle  are  linear  and  independent,  an 
analytical  proof  of  the  validity  of  these  assumptions  would 
require  an  explicit  representation  of  F  and  H. 
Unfortunately,  there  is  no  known  closed  form  solution 
available  to  cover  all  the  possible  sets  of  initial 
conditions.  Therefore,  a  series  of  computer  experiments 
were  performed  to  determine  whether  the  method  of  error 
sensitivity  analysis  was  suitable. 


2 •   Tests  of  the  Proposed  Solution 


a.   Static  Accuracy 

The  first  test  of  the  proposed  solution  was 
designed  to  measure  the  effectiveness  of  the  error  analysis 
technique.  To  separate  the  error  analysis  technique  from 
any  errors  which  might  be  induced  by  airspeed  prediction  or 
dive  anqle  calculation,  it  was  assumed  that  these  functions 
were  performed  exactly  by  the  micro  computer.  To  measure 
the  ability  of  the  system  to  improve  accuracy,  an  arbitrary 
value  of  plus  200  feet  or  minus  200  feet  was  assigned  to  the 
uncorrected  impact  error,  DX  in  figure  6.  The  original  dive 
angle  was  assigned  a  value  which  was  five  degrees  steep  or 
five  degrees  shallow.  It  was  assumed  that  the  correction 
which  was  applied  to  the  gunsight  would  cause  the  bomb  to 
hit  the  target  if  the  pilot  could  establish  a  parallel  dive 
path  whose  displacement  from  the  original  dive  path  was  -DX 
feet  (  measured  along  the  X  axis) .  It  was  assumed  that, 
instead  of  performing  such  a  maneuver,  the   pilot   used   the 
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new  sight  image  to  make  a  single  change  in  his  dive  path. 
This  change  occurred  three  seconds  prior  to  reaching  release 
altitude  and  moved  the  aim  off  point  -DX  feet.  This 
maneuver  causes  a  small  change  in  the  dive  angle  for  which 
the  gunsight  was  not  compensated.  At  release  altitude,  the 
effect  of  this  change  in  dive  angle  was  computed  using  the 
error  sensitivity  tables  for  the  Mk-76  Practice  Bomb  (DX ' 
figure  6) . 

A  FORTRAN  program  was  executed  to  determine  the 
values  of  DX*,  based  on  the  previous  assumptions.  The 
program  is  listed  in  Appendix  C  and  results  are  in  Table  1. 
In  general  the  release  error  sensitivity  method  improved 
accuracy  by  a  factor  of  eight,  and,  in  the  worst  case,  by  a 
factor  of  five. 


31 


CORRECTION  ALTITUDE 


RELEASE  ALTI TUDE 


DIVE   BEFORE  CORRECTING 
DIVE  AFTER  CORRECTING 


.♦.__•      UNCORRECTED   BOMB 

. CORRECTED   BOMB 

Fig    6       SIGHT    ANGLE    CORBECTION     IN    A    DIVE 


32 


r~t  -4  ( \J  cnj  f\)  — « cni  f\i  r\)  r\j  in;  cni  (Xl  cm  cm  cnj  cni  cni  cni  r\  i  cnj  cni'  cm  cni 
OOOOOOOOOOOOOOQOOOOOOOOO 

UJ  LU  LU  LU  UJ  UJ  UJ  UJ  LU  UJ  LU  LU  LU  LLi  Uj  UJ  LU  UJ  LL)  UJ  LU  LLI  LU  LU 

coOh-o— •  4-<j-r-~<7\o>|'Nj-4'oom',oro— "n  a>co;or— -o 
in  ro  -4"  o  o  ro  -4-  ro  n  c/>  m  o  »-h  ro  o  ^  ro  ro  ixi  r-  r—  f— ~h  r- 
r»s)--^roor-fO>i-vtrom-4-0'X)Ocr'0'— io  rocoo-Hr«- 

XCNimsJ-r-cOO  O0  0*0'r»0c0roin^^-r0r"--*<£>v00^ 

«4-  ro  s}-  ro  —*  cr>  c M  i  \J  *-*  — <  cnj  oj  —<  —i  ro  ■  o  ~t  -*  <ni  cnj  cni  cnj  cni  cni 

oooooooooooooooooooooooo 
I      I      I      I      I      I      I      I      i      I      I      I 


—l.-tOOOOOOOOooOOOOOO  — <-HOO— »^ 

oooooooooooooooooooooooo 
I  I 

LU  LU  LU  UJ  LU  IX)  UJ  1 JJ  U  J  LU  'XI  UJ  U  I  LU  UJ  UJ  UJ  LU  UJ  UJ  UJ  LU  LU  LU 

r—  O-4-rOCOrO— <  vJ-iAOOrxi  vOOcomm^OLn  r—  m—trOO 

I— >t  CO  CO  ro  <t"  in  -OOOmp~-cOCNI<r,coroCNjOm  vKNJ—ir-- ^ 

uj  vj-  f-*  <q  — i  o  r-  co  oo  c*  co  o  ^i*-  <r  ro  o  r—  cj>  — *  —*  ro  o  -*j  in 
Xcor\Jr*r\jro— troorxim-JinrNj-j'-tcrNCO^-r-  Q>f-^0>t*« 
H-in— <oo  ^^coDinmo— <tM^OsOooc3or—  _ir»-mcom 
Qr-4rOror<-cMsfvO(\|cooo>>  rocr>or\J^cNjin  — irvi^-OOO 
vO  O  m  in  ro  ro  co  o>  vf  m  co  cp>  ro  >£•  r^  r—  co  co  -*  — <  co  co  — <  •-• 

oooooooooooooooooooooooo 
»      I      I      I      I      I      I      I      I      »      I      I 


2:  co  co  co  co  r\j  r\i  r\j  cnj  m  in  in  ■  r\  in  in  in  in  ■>£>  O  o  -O  -fr  -4-  ^  <*• 
0--i-H~H-^rOrOfOror*-r-r*-  r-^-r-r-r»xjCNJc\iCNJ  coco  coco 

i—  *o  *o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  ro  ro  ro  ro  in  in  in  in 

OC^CT>CJ>Cr>COCOCOCO-t<f*^  vf<f<t- -4"  >!■--*— I— (HNNI^N 

LU<r»rO(^rorMrNir^r\j^w^-<^-i--<--H^(^(^<v0^sOOCT4C7sCr>C)> 
cj:_i,^_ «—«r\jr\jc\jc\isi-sf>a-  -3-m  in  mm  mining  coco  coco 
en 

O  X 

o  o 

Q 

xoooooooooooooooooooooooo  < 

o •• •••••••••• 

oooooooooooooooooooooooo  I— 

oooooooooooooooooooooooo  LU 

r\l  r\)  r\|  cnj  oo  cm  r\j  cni  cm  cnj  r\ic\j  r\l  cnjcnjcnjcnjcnjcnjcni  cm  cm  cm  cm  X 

I  I  I  I  I  I  I  I  f         I  I         I  H- 

uj  o 

> 

i-i  IX. 

o  o 

-J  2 

<oooooooooooooooooooooooo  O 

OOOOOOOOOOOOOOOOOOOOOOOOO  -* 

t— OOOOOOOOOOOOOOOOOOOOOOOO  I— 

OOOOOOOOOOOOOOOOOOOOOOOOO  < 

<oooooooooooooooooooooooo 


inininminininmininminmininmoooommmin 
~^  ft  — <  — *  cnj  cni  cni  cnj  ro  ro  (N)  cnj  ro  ^o  >$■  vj-  in  i  n  in  in  <o  vO 


z> 


00 


> 

— tOOOOOOOOOOOOOOOOOOOOOOOO  -* 
OOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOO  LU 

OOOOOOOOOOOOOOOOOOOOOOOOO  —I 

LUOOOOOOOQOOOOOOOOOOOOOOOO  CO 

a: «••  < 

HOOOoooooooooooooinmmiAOOOO  i— 

LO^^>^,^^^r^rMCNJCNj1-ororOrorO(OrO'nO>^^--4'^--^sOvO^O 

LU 

a 


33 


b.   Eeal  Time  Accuracy 

An  important  part  of  the  proposed  solution  was 
the  calculation  of  the  angle  of  attack  of  the  fuselage 
reference  line.  The  pilot  would  use  this  information  to 
establish  the  aircraft  in  a  straight  path  which  intersected 
the  ground  at  a  precomputed  distance  beyond  the  target, 
called  the  aim  off  point,  reference/"  3  _/.  The  computer 
would  also  use  this  angle  of  attack  information  to  compute 
the  actual  dive  angle  of  the  aircraft. 

Although  all  Naval  carrier  based  jet  aircraft 
are  equipped  with  angle  of  attack  sensors,  these  devices  are 
designed  primarily  to  provide  accurate  information  in  the 
landing  phase  of  flight  and  are  inherently  inaccurate  during 
high  speed  low  angle  of  attack  situations,  reference/"  4  J. 

Therefore,  in  order  to  have  a  more  accurate 
value  for  the  angle  of  attack,  figure  2-56  of  reference/"  2  J 
was  used  as  a  source,  -  Fig  7. 
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Figure  7  is  the  chart  which  a  pilot  would  use  to 
calculate  the  angle  of  attack  for  any  particular  air  to 
ground  delivery  condition.  The  chart  may  be  used  with 
either  of  two  sets  of  input  keys.  The  first  set  is  dive 
angle,  temperature,  true  airspeed,  and  gross  weight.  The 
second  set  consists  of  calibrated  airspeed,  gross  weight, 
and  dive  angle.  Since  the  central  air  data  computer 
provides  calibrated  airspeed  information,  the  second  set  was 
used.  To  calculate  manually  the  angle  of  attack  for  a  given 
calibrated  airspeed,  gross  weight,  and  dive  angle  the 
following  procedure  is  used: 

(1)  Find  the  desired  calibrated  airspeed 
(cas)  on  the  right  half  of  the  bottom  of  the  chart. 

(2)  Move  vertically  to  the  curve  labeled  with 
the  appropriate  weight. 

(3)  Proceed  horizontally  to  the  left  until 
intersecting  the  straight  line  labeled  by  the  desired  dive 
angle. 

(4)  Proceed  directly  down  the  chart  and  read 
the  desired  value  of  the  angle  of  attack  in  mils  from  the 
bottom  of  the  chart. 

To  reproduce  the  information  contained  in  figure 
7,  it  was  necessary  to  represent  the  weight  curves  and 
straight  lines  in  an  algebraic  form.  The  form  chosen  was 
the  polynomial.  Since  all  air-to-ground  charts  are  based  on 
the  assumption  that  the  aircraft  gross  weight  is  42,000 
pounds  with  scalar  adjustments  for  different  weights,  the 
same  scheme  was  adopted  for  the  micro  computer  system.  The 
42,000  pound  curve  was  approximated  with  a  second  degree 
interpolating  polynomial,  F  (cas) .   The  straight   Lines   were 
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considered  to  be  of  the  form 

y=m  *  x  +  b 

where  *  y1  represents  ordinate  values,  'x1  represents 
abscissa  values,  'm'  is  the  slope  of  the  line,  and  'b1  is 
the  *y'  intercept  value.  Substituting  ALPHA,  the  angle  of 
attack,  for  '  x',  ? (cas)  for  •  y',  m (THETA)  for  m,  and 
B  (THETA)  for  'b',  the  following  equation  was  evaluated: 

ALPHA=  (F  (cas) -B  (THETA)  ) /m  (THETA)  . 

F  was  a  quadratic  interpolation  of  the  42,000  pound  curve. 
The  argument  of  F,  cas,  was  taken  as  the  difference  between 
calibrated  airspeed  and  295.  This  linear  translation 
allowed  single  byte  multiplication  to  be  used.  THETA  was 
the  dive  angle.  B  was  a  tabled,  first  divided  difference 
function  cf  the  zero  mil  intercept  values  of  the  straight 
lines  in  the  lower  left  quadrant  of  figure  7.  M,  the  slope 
of  these  straight  lines,  was  found  to  approximate  closely 
the  value  of  '-SECANT  (THETA)  • . 

Therefore,  in  order  to  speed  the  calculations, 
the  equation  which  was  actually  solved  was 

mil  =  cos  (THETA)*  (B  (THETA)  -F  (cas))  .   (13) 

The  COSINE  calculation  was  based  on  a  divided  difference 
table  which  was  arranged  in  four  degree  increments. 

An  ALGOL  program  was  executed  to  simulate  an 
INTEL  8C80  fixed  point  arithmetic  evaluation  of  the 
preceeding  equation.  A  listing  of  the  program  is  contained 
in  Appendix  C.  The  results  of  the  ALGOL  simulation  appear 
in  table  2.  A  ELM  version  of  the  routine  was  timed  by  the 
INTERP8  as  requiring  4798  micro  seconds. 
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c.   Forecasting  Accuracy 

A  third  topic  for  investigation  is  the  ability 
of  the  micro  computer  to  forecast  accurately  the  value  of 
the  release  true  airspeed  early  in  the  dive.  The  basic 
forecasting  scheme  is  to  fit  a  parabola  through  a  series  of 
samples  of  airspeed  and  altitude  .  It  is  assumed  that  both 
values  can  be  read  to  the  equivalent  of  four  decimal  places 
of  precision.  In  fitting  the  parabola,  altitude  is  treated 
as  the  independent  variable  and  the  release  true  airspeed  is 
assumed  to  be 


2 
airspeed  =  C  +  C  *relalt+C  *relalt    (15) 
12  3 

where  'relalt1  is  the  calculated  release  altitude. 


To  determine  whether  such  a  method  was 
appropriate,  it  was  necessary  to  first  construct  a  model  of 
the  F-4  in  a  dive,  and  then  write  a  computer  program  to 
simulate  the  model.  -The  model  chosen  was  that  of  Newton's 
second  law  of  motion: 

F=m*a 
where  force,  F,  was  assumed  to  be 

Thrust  +  W*sin  (THEIA) -drag. 
Thrust  and  drag  were  opposing  forces  acting  parallel  to   the 
flight   path,   THETA   was   constant,   and   thrust,  mass,  and 
weight,  W,  were  also  constant. 

The  resulting  differential  equations  were 

A=  (Thrust  +  K  *  SIN  (THETA)  -  drag) /mass   (16) 
and 
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V=A*t.        (17) 


Drag  was  assumed  to  be  a  function  of  dive  angle 
and  true  airspeed.  Figure  2-61  of  reference  £~  2  J  was 
selected  as  the  source  of  data  from  which  to  obtain  values 
of  drag.  This  particular  figure  represents  true  airspeed 
vs.  altitude  at  various  dive  angles  for  a  particular  engine 
speed,  gross  weight,  and  external  configuration  (drag  index 
60) .  The  relationship  is  represented  by  a  series  of  smooth 
curves,  where  each  curve  corresponds  to  a  particular  dive 
angle.  The  curves  themselves  were  first  represented  using 
an  interpolating  second  degree  polynomial  determined  by  a 
least  sguares  fit  of  eight  evenly  spaced  sets  of  data 
points.  The  coefficients  of  the  polynomial  were  calculated 
by  the  Naval  Postgraduate  School  Computer  Library  subroutine 
LSQPL2.  With  these  coefficients,  it  was  then  possible  to 
determine  fcr  any  altitude  and  dive  angle  a  true  airspeed 
value  based  en  the  curves. 

Using  the  relationship  previously  established, 
it  was  possible  to  determine  an  approximate  value  for  the 
time  reguired  to  fly  from  one  altitude  to  another.  Using 
increments  of  10  knots  and  assuming  that  the  rate  of  change 
of  velocity  was  linear  within  the  interval,  an  empirical 
approximation  for  the  accelleration  in  a  particular  time 
increment  was  calculated.  The  value  of  the  drag  during  that 
time  increment  could  then  be  calculated.  This  method  was 
used  to  establish  pairs  of  data  points  (true  airspeed,  drag) 
for  each  of  the  dive  angles  given.  An  additional  pair  (0,0) 
was  added  for  each  dive  and  the  resulting  data  was  again 
fitted  using  the  routine  LSQPL2.  The  coefficients  were  used 
in  FORTRAN  subroutine  DRAG,  see  Appendix  E. 

With  the  drag  function  computed,  equations   (16) 
and    (17)    were    numerically    solved   using   a   standard 
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Runge-Kutta  fourth  order  integration  scheme.   this  was   done 
in  a  FORTRAN  subroutine,  INTEG,  (see  Appendix  E) . 

Using  the  model,  three  forecasting  schemes  were 
tested.  The  first  scheme  consisted  of  sampling  airspeed  and 
altitude  values  at  one  second  intervals.  A  Newton  divided 
difference  interpolating  polynomial  was  then  generated  and 
release  true  airspeed  was  predicted  as 


tas       =     tas      +    (release    alt-alt   )*(F 

release        t  t      (0.1) 

0  0 

+  (release  alt-alt   ) *F       ).   (18) 

t     (0,1,2) 
1 

F      =(F  -F  )/(x  -x  )  , 
(0,1)    0   1     0   1 

F      is  analagously  defined, 
(1,2) 

F        =  (F      -F      )/(x   -x  )  , 
(0,1,2)     (0,1)   (1,2)     0    2 

X   represents   the    independent    variable    and    F 
represents  the  dependent  variable. 


The  second  method  consisted  of  taking  a  total  of 
twelve  samples  during  a  two  second  period.  The  first  four 
samples  were  added  and  the  average  was  taken  as  the  true 
value.  Similar  averages  were  used  to  represent  the  two 
remaining  inputs  to  a  Newton  divided  difference  second 
degree  interpolating  polynomial  of  the  type  described  above. 

The  third  prediction  scheme  was  to  take  20 
samples  during  the  same  interval  and  to  generate  a  least 
squares  orthogonal  polynomial  to  fit  the  data.  The 
polynomial  was  also  a  parabola  and  it  was  generated  and 
evaluated  based  on  subroutines  contained  in  reference  /5  j. 
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Por  the  experiment,  1000  dive  runs  were 
simulated.   Each  run  was  based  on  the  following  assumptions: 

(1)  Entry  airspeed  equals  345  knots. 

(2)  Thrust  equals  20,000  pounds. 

0 

(3)  Dive  angle  equals  30  . 

(4)  Entry  altitude  equals  6000  feet. 

(5)  True  airspeed  and  altitude  would  be  read  to 
four  decimal  digits  of  precision. 

(6)  On  any  particular  run,  the  altitude  and 
airspeed  error  would  remain  constant.  Between  runs  the 
error  would _  be  normally  distributed  with  the  mean  airspeed 
error  of  zero  and  standard  deviation  of  .5  knot.  The 
altitude  error  would  have  a  mean  of  zero  and  a  standard 
deviation  of  10  feet. 

The  latter  distributions  are  the  result  of 
assuming  that  the  errors  were  normally  distributed  with  the 
maximum  allowable  error  equal  to  three  times  the  standard 
deviation.  Allowable  instrument  errors  were  obtained  from 
ref  erence/6"  y. 

Normal  distributions  were  simulated  using  the 
Naval  Postgraduate  School  Computer  Library  package  LLRANDOM. 
The  test  was  repeated  with  a  different  seed  to  verify 
original  results.   See  Appendicies  A  and  E. 

d.   System  Performance 
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To  determine  whetner  the  use  of  the  micro 
computer  would  significantly  improve  the  delivery  accuracy 
of  the  F-4  in  the  air  to  ground  enviornment,  an  experiment 
was  conducted  based  on  the  following  assumptions: 

(1)  The  previously  described  model  represented 
an  F-4  in  a  dive  bombing  run. 

(2)  The  Mk-8  1  Low  Drag  Bomb  is  a  typical 
ung aided  air  to  ground  weapon. 

(3)  In  the  simulation,  it  would  be  assumed  that 
the  pilot,  at  release  altitude,  would  have  the  sight  aimed 
at  the  target. 

(4)  All  weapons  would  be  released  from  a 
position  7  feet  below  and  18  feet  aft  of  the  gunsight. 


(5)    All   deliveries   would    be    in    no-wind 


conditions. 


(6)  The  path  of  the  Mk-8  1  Low  Drag  Bomb  could 
be  described  by  the  following  set  of  equations  which  were 
obtained  from  reference  //_,/: 


X'  •   =   -E  *  X 


(19) 


and 


i  t 


=   -E  *  Y«   -g   (2  0) 


where  X  represents  down  range  travel,  Y  represents  altitude, 
X* ,  and  Y'  are  the  horizontal  and  vertical  components  of  the 
velocity  vector.  The  value  of  the  other  variables  was 
determined  as  follows: 
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where 


E=  (c  *  p  *  V   *  K  )  /144 
a     d 


c    is    the    reciprocal    ballistic    coefficient 

2 

(diameter  /  weight)  , 

-  (0.00003158  *Y) 
p  is  the  air  density  0.0751366  *  e  , 


V   is  the  true  airspeed  in  ft/sec  .  , 
a 


and 


K   is  the  drag  coefficient 
d 


The  experiment  consisted  of  simulating  1000  dive 
runs  for  each  of  the  following  standard  deliveries: 

0 
20   dive,  2000  ft.  release,  450  knots  true  airspeed 

0 
30   dive,  3000  ft.  release,  450  knot  true  airspeed 

0 
30   dive,  5000  ft.  release,  500  knot  true  airspeed 

0 
45   dive,  4000  ft.  release,  450  knot  true  airspeed 

0 
45   dive,  6000  ft.  release,  500  knot  true  airspeed 

0 
60   dive,  7000  ft.  release,  500  knot  true  airspeed. 


On  each  run  the  sampling  method  previously 
described  as  one  second  interval  sampling  was  used.  The 
predicted  airspeed  was  calculated  by  equation  (18).  The 
entry   conditions  were  considered  to  be  normally  distributed 
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about  a  correct  set  of  (airspeed,  altitude,  thrust,  and  dive 

angle)   which  was  calculated  to  achieve  the  desired  delivery 

parameters.   The  standard  deviations  of  the  distributions  of 

the   entry   conditions   were  10  knots  for  true  airspeed,  200 

0 
feet  for  altitude,  2.5    for  dive  angle,  and  750  pounds   for 

thrust. 

Normal  distributions  were  simulated  using  the 
Naval  Postgraduate  School  subroutine  LLRANDOM.  The  standard 
deviations  were  assigned  with  a  data  statement.  Two 
pseudorandom  variables,  x  and  y,  were  obtained  from  the 
LLRANDOM  package.  The  variable  x  was  obtained  from  a 
uniform  (0,1)  distribution  and  y  was  obtained  from  a  normal 
(0,1)  distribution.  The  desired  value,  z,  was  then 
calculated  as 

z=  DMEAN  +  y  *  STDV  *  COIN (x)  where  DMEAN  is  the 
desired  mean,  STDV  is  the  standard  deviation,  and 
COIN  is  a  statement  function  which  returns  a  +1  if 
x<.5  and  a  -1  otherwise. 

After  the  release  conditions  had  been  predicted 
by  the  subroutine  SAMP,  two  conditions  were  simulated. 
First,  that  the  aircraft  had  no  micro  computer  and  that  the 
pilot  continued  in  a  steady  dive  to  release  altitude. 
Second,  that  the  pilot  had  a  micro  computer,  that  the  micro 
computer  supplied  his  gunsight  with  a  corrected  sight 
setting  egual  to  the  value  calculated  by  SAMP  ,  and  that  the 
pilot  adjusted  his  aircraft's  dive  path  as  soon  as  this 
information  was  available  ,  so  that  at  release  altitude,  the 
new  sight  setting  was  'on'  the  target. 

For  each  condition,  a  simulated  bomb  was  dropped 
at  release  altitude.  A  miss  for  the  bomb  was  calculated  as 
the   distance   between   the   point   where   the   sight   angle 
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intersected  the  earth  and  the  down  range  travel,  x,  which 
was  numerically  computed  (  based  on  the  equations  (19)  and 
(20)).  In  addition,  in  the  case  of  the  modified  flight 
path,  a  correct  sight  setting  was  also  calculated.  Finally, 
the  ratic  of  the  miss  distance  with  the  micro  computer  to 
the  miss  distance  without  the  micro  computer  was  calculated. 
If  the  micro  computer  gave  a  greater  miss  than  the  basic 
aircraft,  the  ratio  was  set  to  plus  or  minus  1.1,  depending 
on  whether  the  errors  were  both  on  the  same  side  of  the 
target. 

Arrays  were  built  with  these  1000  values  and 
histograms  were  printed  showing  the  distribution  of  the 
various  values.  See  appendix  B.  The  FORTRAN  routine  which 
simulated  the  experiment  is  in  appendix  F.  An  additional 
table  was  constructed  from  the  histograms  which  compares  the 
median  miss  distance  under  the  assumptions  that  the  micro 
computer  *as  or  was  not  used.   See  table  3. 
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e.   Micro  Computer  Suitability 

The  previously  tested  error  sensitivity 
algorithm  was  converted  from  FORTRAN  to  PLM  for  test  on  the 
INTEL  8080-  The  program  was  compiled  on  the  Naval 
Postgraduate  School  I3M  360/65  and  interpreted  using  the 
INTSHP8  subroutine.  Timing  data  was  obtained  and  found  to 
be  .75  seconds  for  the  entire  program.  See  Appendix  G.  The 
program  was  then  executed  on  the  INTEL  8080  and  final  output 
(sight  setting)  was  found  to  be  in  agreement  with  sample 
results  which  had  been  obtained  from  executing  the  FORTRAN 
version  on  an  SDS  9300  computer. 


B.    AIR  TO  AIR 


1 .   General 

In  air-to-air  gunnery  considerations,  it  must  be 
assumed  rhat  the  position  of  the  target  is  some  function  of 
time, 

(x,y,z)  =  F  (a{t)  ,b<t)  ,c(t))  . 

Since,  in  general,  no  assumptions  can  be  made  to 
restrict  this  motion  relative  to  the  (x,y,z)  coordinate 
system  of  air  to  ground  considerations,  a  new  coordinate 
system  is  introduced.  The  origin  of  the  new  system, 
{x1  ,y*  ,2.')  is  the  muzzle  of  the  fighter  aircraft's  cannon. 
The  x'  axis  is  always  oriented  so  that  a  bullet's  initial 
velocity  vector  will  be  aligned  with  the  positive  x'  axis. 
The  y'  axis  is  perpendicular  to  the  x'  axis  and  forms,   with 
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the  x'  axis,  a  plane  which  contains  the  wing  roots  of  the 
fighter.  The  z'  axis  is  orthogonal  to  the  (x',y')  plane 
with  the  positive  direction  measured  in  the  same  direction 
as  the  vertical  tail  of  the  fighter. 

It  should  be  noted  that  the  x1  axis  is  not  always 
aligned  with  the  gun  barrel  boresight  line.  This  is  caused 
by  the  bullet's  velocity  vector  being  the  vector  sum  of  the 
basic  muzzle  velocity  vector  and  the  aircraft  velocity 
vector.  When  the  aircraft  is  flying  at  an  angle  of  attack 
greater  than  zero  these  two  vectors  are  not  coincident. 
Therefore,  the  x'  axis  is  shifted  away  from  the  gun 
boresight  line  toward  the  relative  wind. 

The  set  of  points  from  which  a  fighter  can 
effectively  employ  guns  against  an  airborne  target  lies 
within  a  30  degree  cone  emanating  from  the  tail  of  the 
target.  The  minimum  firing  range  within  this  cone  is 
approximately  500  feet  and  the  maximum  range  is 
approximately  2000  feet. 

The  dimensions  of  this  firing  region  effectively 
restrict  the  problem  so  that,  for  a  bullet  to  hit  the 
target,  the  time  of  flight  is  less  than  two  seconds.  This 
short  time  of  flight,  together  with  the  basic  muzzle 
velocity  of  the  standard  aircraft  cannon  (3300  fps) ,  allow 
the  computations  of  the  proper  lead  angle  to  be  performed  in 
four  distinct  steps,  Reference  /8  J7.  The  steps  are  as 
follows: 

(1)  Estimate  the  future  angular  displacement  of 
the  target  from  the  present  line  of  sight.  This 
displacement  is  assumed  to  be  caused  by  the  target's  present 
velocity  vector  acting  as  a  constant  during  the  time  of 
flight  of  the  bullet.  This  angle  is  called  the  lead  for 
target  motion. 
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(2)   Calculate  the  effect  of  gravity  drop. 


(3)   Calculate  the  effect  of  the  fighter's  angle 


of  attack. 


(4)  Calculate  the  effect  of  the  distance 
between  the  gun  and  the  gunsight  (parallax) . 

2 .   Lea d  for  Target  Motion 

The  calculation  of  the  lead  for  target  motion  in  a 
plane  parallel  to  the  (y'/Z1)  plane  is  based  on  the 
following  analysis  of  the  situation: 

(1)  If  the  fighter  is  tracking  the  target, 
that  is,  the  fighter  is  in  the  firing  cone  with  the  35  mil 
pipper  held  steady  on  the  target,  the  relative  positions  of 
the  two  aircraft  are  'close'  to  the  positions  needed  for 
successful  gun  firing. 

(2)  Therefore,  assume  that  the  sight  line  and 
the  x*  axis  are  coincident. 

(3)  At  this  time,  the  fighter's  velocity  minus 
the  rate  of  closure  represents  the  target's  velocity 
component  parallel  to  the  x'  axis. 

(k)  The  rate  at  which  the  fighter  is  turning, 
in  radians/second,  multiplied  by  the  range  is  a  close 
approximation  to  the  target's  velocity  component  parallel  to 
the  (y ■ , z ')  plane. 

(5)  If  there  exists  a  time,  t',  such  that  a 
bullet   fired  at  time  t=0  will  strike  the  target  at  time  t*, 
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then  the  geometry  of  the  flight,  of  the  bullet  in  a  plane 
containing  the  x'  axis  and  the  component  of  the  velocity 
vector  of  the  target  which  is  in  the  plane  parallel  to  the 
(y'/Z')  plane,  can  be  approximated  by  Fig  8. 
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Fig  8    AIR  TO  AIR  GUNNERY  DIAGRAM 
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(6)  Therefore,  the  lead  for  target  motion 
orthogonal  to  the  x'  axis  can  be  approximated  by  the  angle, 
psi,  (see   Fig  8) . 


The   calculation   of  the  angle  in  mils  uses  the  fact 


that 


lim    sin  (x) /x  =1.   This  implies  that,  for  small  angles, 
x-0 


ARCSIK  (x/y)  =   x/y  = 


AECTAN  (x/y) 


where  =  means  "is  approximately  equal  to". 

In  fact,  for  angles  smaller  than  200  mils,  the  maximum  error 
of  these  approximations  is  4  mils.  Since  the  set  of 
possible  answers  lies  within  the  range  0-200  mils,  the 
following  equation  was  solved  to  determine  psi. 

psi=  VO/VX1. 
VQ,  the  component  of  the  target's  velocity  vector  which  was 
orthogonal  to  the  X'  axis,  was  calculated  as  the  product  of 
the  fighter's  rate  of  turn  multiplied  by  the  range  to  the 
target,  multiplied  by  1000.  VX',  the  average  velocity  of 
the  bullet,  was  calculated  as  the  range  to  impact  divided  by 
the  time  of  flight. 

Since  the  calculation,  at  this  stage,  ignores 
gravity,  and  since  the  bullet  otherwise  satisfies  the 
assumptions  of  the  air  to  ground  weapon  ,  provided  that  the 
x'  axis  is  parallel  to  the  surface  of  the  earth,  equation 
(19)  can  be  modified  to  give 


=   -E  *  V 


(21) 
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Since   the  velocity  of  the  bullet  during  it's  time  of 

flight  will  be  several  times  larger  than  the  speed  of  sound, 

and   since  the  shape  of  the  K   vs.  Kach  curve  is  almost  flat 

d 

above  1.6  times  the  speed  of  sound,   equation   (21)   can   be 
further  simplified  to 


-p*  K  *V 


(22) 


or  equivalentiy,  since  V    >  0, 


/  V 


=  -p  *  K 


Since  the  effect  of  gravity  is  ignored,  the  density,  p,  will 
remain  constant  so  that  (22)  becomes 


/  V 


2 


=  -K  ,  a  constant.  (23) 


Integrating  both  sides  with  respect  to  t  gives 

-1/V     =  -K   *  t   +   C  .   (24) 
x'       2  1 

Solving  (24)  for  the  initial  conditions  of  the  bullet  as   it 

leaves  the  gun  gives 


1/V 


(K   *  V   *  t   +1  ) /V    (25)  where  V   was   taken 
2     m  m  m 


to   be   3300  *  V  ,  and  V   was  the  fighter's  velocity  in  feet 
per  second.   Expressing  (25)  in  an  equivalent  form  yields 


V/(K   *  V   *  t   +   1).  (26) 
m    2     m 


Integrating  equation  (26)  again  with   respect   to   t 
gives  the  equation  for  distance  traveled  by  the  bullet 
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X'   =   ln(K   *  V   *  t   +   1)/K    +   C  -     (27) 

2     m  2       3 

Since  x'  at  t=0  is  zero,  C  =0. 


The  time  of  flight  may   now   be   calculated   as   the 

time,  t    ,  when 
tof 


RANGE   +   {V    -   V  )  *  t 

f       c      tof 

ln(K   *  V   *  t      +   1)/K  .   (28) 
2     m     tof         2 

RANGE  is  the  distance  to  the  target  at  time  t=0.   V   is   the 

c 

closing  velocity  at  time  t=0. 

Since  the  argument  of  the  natural  logarithm  in 
equation  (28)  is  strictly  greater  than  1,  a  Newton's  fixed 
point  iteration  function  may  be  used  to  find   the   value   of 

t    .   The  iteration  function  is 
tof 


t     =  t   -  F(t  )/F«  (t  )  .   (29) 
n  + 1     n      n      n 

F(t)  and  F1  (t)  were  defined  as  follows: 


F(t)  =  RANGE  +  (V   -  V  )  *  t   -In  (V   *  K   *  t  +  1) /K    (30) 

f     c  m     2  2 


and 


F«  (t)  =  (V   -  V  )   -  V  /(V   *  K    *  t   +   1).   (31) 

f    c      m   m     2 


The  first  approximation  of  t     was  used  to  find  the 

tof 

values  VX*,   and   subsequently   psi.    To   reduce   the   time 
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required   to   calculate  t  ,  t   was  always  assigned  the  value 

2    1 


one . 


3-   The  Effect  of  Gravity  Dro£ 

In  calculating   the   lead   for   target   motion,   the 

effect   of  gravity  was  ignored.   A  true  trajectory,  however, 

would  always  curve  toward  the  earth  as  described  by  equation 

(20).    For  a  time  of  flight  less  than  two  seconds,  equation 

<20)  may  be  simplified  by  assuming  that  E  =  0.   This  allows   a 

closed   form   solution   for   the   distance   travelled  by  the 

bullet  below  the  no  gravity  flight  path  during   time/   t    , 

tof 

seconds.   The  solution  is 


2 
y   =   32.  17  *  t  /2      (32) 

The  error  in  (32)  is  less  than  20  feet  reference  f  8  J '. 


The   angle,   in   milliradians,   subtended    by    the 

distance  y    ,  at  a  future  range  of  x1     feet  is 
tof  tof 


mils  =y     *  1000/x1 

gravity  drop      tof  tof 


**  •   Effect  of  Anqle  of  Attack 

As  the  bullet  leaves  the  barrel  of  the  gun,  it 
aligns  with  the  relative  wind.  That  is,  the  velocity  vector 
of  the  bullet  shifts  from  the  gun  boresight  line  to  the  x' 
axis.   The  amount  of  this  shift,  in  mils,  is  calculated  as 
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MILS     =  V   *  ACA/(V   +  3300)  , 
AOA     f  f 


where   AOA   represents   angle  of  attack,  in  mils,  of  the  gun 
barrel  boresight  line. 


5.   Parallax  Correction 


Since  the  gun  and  the  sight  pivot  are  approximately 
seven  feet  apart  in  the  z*  direction,  a  small  factor  must  be 
added  to  the  sight  setting.  In  the  F-4,  this  setting  is 
3.19  mils,  reference  /Q/  .  Fig  9  shows  that  the  effect  of 
adding  a  parallax  correctio-n  is  to  have  the  zero  sight  line 
(ZSL)  intersect  the  gun  boresight  line  (TBL)  at  a  known 
distance  (harmonization  range)  in  front  of  the  fighter. 


GU*  80RES/GH 


r 


Fig  9    PARALLAX  CORRECTION 

The  calculations  described  above  were  coded  into  a 
PLM  subroutine  and  tested  on  the  INTERP8  simulation 
subroutine.  The  time  of  execution  was  found  to  be  less  than 
.2  seconds.  In  addition,  a  simulated  sampling  and  data 
preparation  routine  was  timed  at  approximately  ,15  seconds. 
See  appendix  H. 
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CONCLUSIONS 


Based  on  the  experiments  performed  and  the  timing  data 
obtained,  the  micro  computer  is  functionally  capable  of 
providing  an  improvement  in  air  to  ground  delivery  accuracy, 
which  is  on  the  order  of  a  two  fold  increase  in  accuracy. 
Since  the  cost  of  micro  computers  in  decreasing,  whereas  the 
cost  of  fuel  and  replacement  aircraft  is  increasing,  and 
since  improved  delivery  accuracy  means  fewer  sorties  and 
aircraft  exposed  tp  hostile  fire,  a  program  to  study  the 
engineering  factors  and  actual  airborne  effectiveness  of  the 
proposed  micro  computer  modification  is  warranted. 

The  air  to  air  gunnery  implementation  was  not  compared 
to  the  fixed  sight  solution  since  there  is  no  criteria  by 
which  to  judge  the  accuracy  of  the  pilot's  guess  in  the 
latter  case.  The  reaction  time  of  the  computer  (less  than 
.5  second)  appears  to  be  adequate. 

Althouqh  no  live  tests  were  conducted  on  the  system  in 
either  mode,  such  tests  could  be  conducted  at  a  relatively 
low' cost  by  programming  an  A-6E  or  A-7E  airborne  computer  to 
simulate  the  accuracy  and  timing  of  the  proposed  micro 
computer  system.  This  would  eliminate  the  need  co 
reconfigure  the  hardware  of  the  test  aircraft,  at  least  in 
the  preliminary  stages. 

The  cost  of  the  micro  computer  and  associated  10,006 
bytes  of  memory  necessary  to  execute  the  air  to  ground  and 
air  to  air  programs  is  one  thousand  one  hundred  and  fifty 
dollars   based   on  currently  available  prices.   Overall,  the 
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micro  computer  represents  one  way  to   help   curb   the   rapid 
increase  in  the  cost  of  aircraft  procurement. 
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APPENDIX    C 


STATIC    ACCURACY    SIMULATION    PROGRAM 
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APPENDIX  D 


SIMULATION  OF  ANGLE  OF  ATTACK  MODULE 
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APPENDIX  F 
DELIVERY  ACCURACY  COMPARISON 
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